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Purpose: Primary angle-closure glaucoma (PACG) in dogs is usually caused by the gradual collapse of the iridocorneal 
angle and cleft, eventually leading to aqueous humor (AH) outflow obstruction. The condition occurs in several breeds of 
dogs and the prognosis for affected animals is typically poor. We have identified several basset hound (BH) pedigrees, as 
well as unrelated cases with characteristic PACG that in many aspects recapitulates PACG in human patients. The goal 
of this study was to utilize the BH PACG model to characterize the genetics of PACG, and potentially discover genetic 
factors contributing to PACG in humans and animals. 

Methods: We conducted a genome-wide logistic regression test for association using 37 PACG cases and 41 unaffected 
controls. Population stratification and cryptic relatedness were assessed using a multidimensional scaling analysis. 
The expression of two candidate genes within the target tissues of the BH eye was assessed by immunohistochemistry. 
Results: We report significant associations at two novel loci, specifically BICF2P31912 in COL1A2 on chromosome 14 
with a per-allele odds ratio (OR, 95% confidence interval [CI]) of 3.35 (1.73-6.51), P geiiome =3.6xl(r 4 ; and BICF2P893476 
residing in proximity to RAB22A on chromosome 24 with a per-allele OR (95% CI) of 3.93 (1.78-8.66), ^ genome =4.9xl0" 4 . 
COL1A2 and RAB22A demonstrated widespread expression throughout the eye and were prominently noted in the 
ciliary body (CB), trabecular meshwork (TM), and iris. 

Conclusions: Our finding of two genetic associations supports the potential segregation of PACG risk-conferring vari- 
ants in the BH. The genetic associations identified may contribute to mechanisms underlying the pathogenesis of PACG, 
which remain to be elucidated. 



Glaucoma is the third most frequent cause of blind- 
ness worldwide [1], and is known to affect both humans 
and numerous dog breeds. Several types of glaucoma exist, 
including primary angle-closure glaucoma (PACG), which is 
three times more prevalent in humans of Asian descent than 
in European populations [2-4]. In human patients, PACG is 
typically caused by the collapse of the iridocorneal angle due 
to the anterior movement of the iris root, which in turn blocks 
the outflow of aqueous humor (AH) through the trabecular 
meshwork (TM). Conversely, clinical investigation of heredi- 
tary PACG in the basset hound (BH) has revealed progressive 
collapse of the ciliary cleft (CC) and TM region in association 
with a gradual increase in intraocular pressure (IOP), which 
becomes statistically significant (20-24 mmHg) by 20 months 
of age [5]. The condition is marked by retinal ganglion cell 
death, retinal neuroinflammation, axonal degeneration, and 
cupping of the optic nerve head secondary to increased IOP 
[6]. 

BHs display glaucomatous neuropathy similar to that 
observed in human PACG patients. Progressive deficits in 
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retinal ganglion cell function are noted by pattern electroreti- 
nography (pERG) in the basset. Additionally, axonal degen- 
eration and prominent structural changes to the optic nerve 
head are observed in the advanced stages of the disease. The 
frequent absence of early symptoms makes PACG difficult to 
detect and places a large proportion of both human and canine 
cases at risk for blindness. 

There is strong evidence that human PACG is a complex 
disease with a significant underlying genetic component 
[7,8]. Family history is one of the major risk factors for PACG 
and heritability estimates have shown a 3. 7- fold increase in 
disease risk among siblings compared to the general popu- 
lation [9,10]. Despite extensive research efforts, the genetic 
basis of PACG is still incompletely understood. Several 
candidate genes for human PACG have been identified so 
far. A recent genome-wide association study (GWAS) has 
reported three susceptibility loci in an Asian population [8]. 
Association studies have also suggested the involvement 
of several single-nucleotide polymorphisms (SNPs) in the 
matrix metalloproteinase-9 (MMP9; OMIM 120361) in both 
Chinese and Caucasian populations [11,12]. Some of these 
findings were not reproducible in other ethnic populations, 
however, which suggests the genetic heterogeneity of PACG 
in human populations [13]. 
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Similarly, data regarding the genetics of glaucoma in 
dogs is sparse. Genetic investigation of primary open-angle 
glaucoma (POAG) in a beagle colony revealed a single locus 
on chromosome 20 harboring the gene ADAMTS10 [14]. A 
GWA analysis of a late-onset form of PACG described in a 
Dandie Dinmont terrier cohort has also led to the identifica- 
tion of a novel susceptibility locus on canine chromosome 8 
[15]. The locus shares synteny to a region on human chromo- 
some 14q, which harbors several genes associated with POAG 
and primary congenital glaucoma [16]. 

BHs are among several breeds of dogs that are predis- 
posed to developing glaucoma [17]. A total of 5.44% of all 
BHs treated in teaching hospitals presented with glaucoma, 
representing a significantly higher fraction than most breeds 
[18]. Unfortunately, treatment options are limited and affected 
animals often develop bilateral blindness. Other breeds likely 
to be affected include the American cocker spaniel, wire fox 
terrier, and Boston terrier [18]. The observation of multiple 
breeds at risk for developing the same form of the disease 
may possibly indicate that a shared founder mutation has been 
preserved with selective breeding. If true, the identification 
of a significant PACG risk-conferring variant in the BH may 
serve to determine the risk of disease in other breeds sharing 
the same mutation. 

We have identified several BH pedigrees, as well as unre- 
lated dogs with characteristic PACG that closely recapitulates 
the phenotype observed in human patients. All affected dogs 
present with a comparable disease phenotype. Our goal was 
to identify valid genetic associations that confer risk to PACG 
by using the BH as a genetic model. The findings of this 
study are anticipated to provide valuable insight into the 
pathophysiology and genetic mechanisms underlying canine 
PACG. The identification of genetic risk variants in dogs 
can assist in reducing the incidence of PACG in susceptible 
breeds using selective breeding. 

METHODS 

Animals and clinical presentation: All animal studies were 
conducted in accordance with the Association for Research 
in Vision and Ophthalmology (ARVO) statement for use of 
animals in ophthalmic and vision research. The procedures 
conducted were approved by the respective University of Iowa 
and Iowa State University Committees on animal care and 
use. Our GWA investigation of PACG in cases versus controls 
included unrelated as well as related animals derived from 
six pedigrees with clinically confirmed PACG. All animals 
included were examined by a veterinary ophthalmologist as 
described previously [5]. 
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Examination of the anterior segment was performed 
by gonioscopy using a gonioscopy lens (Koeppe; Inc.), 
a hand-held Kowa SL-15 slit lamp (Kowa Optimed Inc., 
Torrance, CA) and a high-frequency ultrasonography unit 
with a 35-MHz probe (E -Technologies, Bettendorf, IA). 
Additionally, gonioscopy photography was performed using 
a gonioscopy lens and a handheld retinal camera (RetCam; 
Massie Research Laboratories, Pleasanton, CA). Collapse 
of the iridocorneal angle was considered the most signifi- 
cant determinant of glaucoma diagnosis. Examination of 
the anterior segment was performed by gonioscopy using a 
gonioscopy lens (Koeppe; Ocular Instruments Inc., Bellevue, 
WA), a hand-held Kowa SL-15 slit lamp (Kowa Optimed 
Inc., Torrance, CA) and a high-frequency ultrasonography 
unit with a 35 -MHz probe (E-Technologies, Bettendorf, IA). 
Additionally, gonioscopy photography was performed using 
a gonioscopy lens and a handheld retinal camera (RetCam; 
Massie Research Laboratories, Pleasanton, CA). Animals 
with possible secondary glaucoma or other confounding 
ocular conditions were excluded from the study. Control 
animals were at least four years of age, and had a normal 
iridocorneal angle and cleft conformation observed by goni- 
oscopy and high frequency ultrasound to ensure accurate 
diagnosis of their unaffected status. Additionally, globes 
(n=6) from animals that underwent enucleation due to clinical 
glaucoma refractive to medical or surgical treatment were 
submitted for histopathology and further evaluation. 

Genetic analysis: Whole blood was collected from 37 clini- 
cally confirmed glaucoma cases and 41 unaffected controls 
with the owners' consent. Genomic DNA was extracted from 
EDTA-stabilized whole blood samples using the Qiagen 
DNeasy blood and tissue kit (Qiagen, Hilden, Germany). 
DNA samples were eluted with deionized water and stored 
at -20°C. Seventy-eight samples were genotyped using the 
Illumina CanineHD BeadChip (Illumina, San Diego, CA), 
which contains 172,000 markers placed on a CanFam2.0 
reference sequence. 

Prior to GWA analysis, preliminary quality control 
assessment of genotyped SNPs was conducted using PLINK 
[19] to remove genotype errors and uninformative data. All 
SNPs were subjected to strict quality control (QC) criteria: 
SNPs with a minor allele frequency <1%, a call rate <95%, a 
rate of missing genotype >10%, and those that departed from 
Hardy-Weinberg equilibrium (P value <0.001) were excluded 
from the analysis. All samples had less than a 5% missing 
genotype rate. 

Pairwise clustering based on the GWA proportion of 
alleles shared identical by state (IBS) between any two indi- 
viduals was performed to account for increased relatedness. 
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All sample pairs were ensured to display an IBS relatedness 
value (pi-hat) <0.9. Population stratification was assessed 
using a multidimensional scaling analysis (MDS), with three- 
dimensional components extracted using PLINK [19]. MDS 
cluster analysis included a selected number of individual 
animals from each pedigree. To our knowledge, our pedigrees 
do not share a recent common ancestor and have not been 
recently interbred. A logistic regression test for GWA was 
performed. GWA testing was conducted with the inclusion 
of MDS-based covariates to correct for increased relatedness 
and population stratification by clustering individual samples 
based on their pedigrees of origin. The genomic-control 
inflation factor (kgc) was applied to the test statistic following 
MDS adjustment by three components. The odds ratio (OR) 
and 95% confidence interval (CI) were estimated for an addi- 
tive effect model of allele dosage in the context of logistic 
regression. The threshold for genome -wide significance was 
set by permutation testing using 1,000 permutations. For 
haplotype definitions we performed linkage disequilibrium 
(LD)-clumping (settings; r2=0.8, pl=0.001, p2=0.01, distance - 
kb=500) using PLINK vl.07 [19]. Quantile-quantile (Q-Q) 
and Manhattan plots of-log 10 p values were constructed using 
RStudio (2012). 

Tissue processing and histological analyses: Eyes derived 
from affected and unaffected animals were fixed immediately 
following enucleation in 4% paraformaldehyde. Tissues were 
dehydrated then subjected to paraffin embedding using stan- 
dard procedures. Globes were embedded and 4 um sagittal 
sections were made using a rotary microtome. Tissue sections 
were deparaffinized in xylene for 20 min then hydrated in 
decreasing ethanol concentration for 5 min at each concen- 
tration. Heat-assisted antigen retrieval was performed using 
lOmM citrate buffer (pH 6.0) for 10 min at 95-100 °C. 
Blocking of endogenous peroxidase activity was achieved by 
treating sections with 2% H 2 0, -methanol solution for 20 min 
at room temperature. Sections were blocked in 5% normal 
goat serum/1% bovine serum albumin (BSA) solution in IX 
Tris-buffered saline (TBS) for 1 h at room temperature then 
incubated overnight at 4 °C in primary rabbit anti-COLlA2 
(1:100) or rabbit anti-RAB22A (1:100) in IX TBS (Protein- 
tech™, Chicago, IL). Sections were treated with biotinylated 
secondary anti-rabbit immunoglobulin G (IgG) antibody 
(1:100; Vector, Burlingame, CA). Color development was 
achieved using an Avidin/Biotinylated-peroxidase Complex 
VECTASTAIN® ABC kit (Vector, Burlingame, CA) by 
direct application to tissue sections for 2 min. Sections were 
dehydrated, cleared in xylene and mounted using Permount™ 
mounting media (Thermo Fisher Scientific Inc., Waltham, 
MA). Immunohistochemical staining was evaluated at the 
light microscopic level. 
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To evaluate collagen localization, tissue sections 
were stained with Picrosirius red stain (Polysciences Inc. 
Warrington, PA) and counterstained with Mayer's Haema- 
toxylin stain (Sigma-Aldrich®, St. Louis, MO) using standard 
procedures. Negative control tissue sections were generated 
by staining with Mayer's Haematoxylin stain only. Picrosirius 
red-stained sections were examined using polarizing micros- 
copy and a color wavelength (from yellow-orange to deep red) 
indicative of collagen fiber thickness was noted in collagen 
containing regions. 

RESULTS 

Animals and clinical presentation: We previously described 
that the development of elevated IOP and subsequent vision 
loss is correlated to the gradual narrowing of the irido- 
corneal angle, resulting in complete angle closure in all 
affected animals. In concordance with these findings, all 
cases included in this study displayed completely collapsed 
iridocorneal angles and CCs confirmed using high-resolution 
ultrasound (HRUS) and gonioscopy examination. Angle 
collapse was observed in association with notable IOP eleva- 
tion of at least 25 mmHG All unaffected control animals 
displayed normal IOPs, and had normal appearing iridocor- 
neal angles and CCs by gonioscopy and HRUS, respectively. 

Genome-wide association mapping: We performed GWA 
analysis of PACG in the BH by analyzing clinically confirmed 
cases and controls to identify genetic variants associated with 
PACG. The genotyping success rate was set to a value >90%. 
Prior to GWA analysis, a strict QC criterion was applied to all 
genotyped SNPs. Out of 172,105 SNPs, we excluded 16,073 
(9.33%) uncalled genotypes, 201 (0.12%) markers due to the 
departure from Hardy-Weinberg equilibrium, 5,739 (3.33%) 
due to a call rate below 95%, and 37,636 (21.87%) non-infor- 
mative markers due to a minor allele frequency below 1%. In 
total, 112,456 SNPs were included in the analysis following 
QC testing. 

All samples analyzed for the presence of excess propor- 
tion of alleles shared IBS revealed a pi-hat value <0.9 and 
were therefore included in the analysis. To account for 
population stratification and reduce the confounding effect 
of subject relatedness, an MDS analysis of all 78 samples was 
conducted. The utilization of pedigree-based MDS clustering 
of samples as a test statistic covariate has proven successful 
at reducing the confounding effect of relatedness in other 
research studies and as such was adopted in our analysis 
[20,21]. A scatter plot revealing pedigree-based clusters was 
generated for 3 MDS components (Figure 1). The MDS plot 
illustrates the formation of six clusters, which appeared to be 
influenced by the samples' pedigrees of origin, as follows: 
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Figure 1. Multidimensional scaling 
analysis (MDS) scatter plot of cases 
and controls plotted for the first 
three MDS dimensions to correct 
for relatedness and population 
stratification. Each dot represents 
a specific dog. Related and unre- 
lated samples are shown clustered 
according to their affiliation with 6 
basset hound (BH) primary angle- 
closure glaucoma (PACG) pedi- 
grees. The group of unrelated dogs 
is seen randomly distributed among 
all pedigree derived samples. 



pedigree 1: 21 samples, 14 unaffected, 7 affected; pedigree 2: 
5 samples, all affected; pedigree 3: 2 samples, all unaffected;, 
pedigree 4: 5 samples, 2 unaffected, 3 affected; pedigree 5: 
3 samples, 2 affected, 1 unaffected; pedigree 6: 17 samples, 
8 unaffected, 9 affected; and unrelated samples: 25 samples, 
13 unaffected, 12 affected. Overall, our samples comprised 
28 males (14 unaffected, 14 affected) and 50 females (27 
unaffected, 23 affected). The group of unrelated samples 
appeared randomly distributed among the six-pedigree- 
derived sample clusters. Logistic regression GWA testing was 
conducted with the inclusion of MDS-generated covariates. 
An inflation factor (X =1.30, y 2 =1. 33) was calculated 

v no correction ' ' 

before covariate adjustment. Q-Q scatter plots of the -loglO 
(p values) expected under the null hypothesis of no genetic 
association versus the observed -loglO (p values) are shown 
for all subjects before (A) as well as following correction 
for sample relatedness (B; Figure 2A,B). In the absence of 
covariate correction, significant deviation from normal 
distribution was observed for the generated p values, which 
demonstrates the strong effects of population stratification 
that is commonly observed when sampling from a founder 
population (Figure 2A). Conversely, covariate adjustment by 
the third MDS dimensional component resulted in a corrected 
inflation factor (\ mectei =h f =0.93). The scatter Q-Q plot 
in Figure 2B reveals significant improvement in the initially 
observed deviation from normal distribution of all p values 



generated following covariate correction and illustrates the 
validity of the method used here to correct for relatedness and 
population stratification. 

We identified two regions in association with PACG 
in the BH following multiple hypothesis testing correction 
of raw p values by permutation. The first and most statis- 
tically significant association (P =0.00036, OR=3.35, 

J ° v genome 

CI 95% =1.73-6.51) was identified for marker BICF2P31912 
on chromosome 14, residing in an intronic region of 
COL1A2 (Figure 3). The marker is part of a four SNP locus 
(Chrl4:19,911,001-20,161,348bp) spanning 0.25Mb. The 
four markers constituting the haplotype displayed complete 
linkage disequilibrium (LD; 23—1, r 2 =l). Notable difference in 
allele and genotype frequency distribution for BICF2P31912 
was observed among cases and controls, (58.1% "G") in 
cases versus (29.3% "G") in controls (Table 1). The allele and 
genotype distribution of BICF2P31912 is representative of all 
SNPs within the locus identified. 

To define the associated haplotypes, we performed an 
LD-based clumping analysis with an LD value of r 2 =0.8. 
Using a threshold value of significance (0.0001 < p<0.001), 
the four previously identified markers of BICF2P31912, 
BICF2S2309101, TIGRP2P188010_rs8723846, and 
BICF2P3 16805 were validated for constituting an LD-haplo- 
block while reaching the highest P values relative to all 

° ° genome 
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Figure 2. Quantile-quantile (Q-Q) 
plot of primary angle-closure 
glaucoma (PACG) genome-wide 
association (GWA) analysis. A: 
Scatter plot of the -loglO (p values) 
expected under the null hypothesis 
of no genetic association versus 



A ExpMed (-109P) B Expected (*gp) the observed -loglO (p values) for 

all subjects before correction for 

population stratification using multidimensional scaling analysis (MDS) generated covariates, and (B) following correction for population 
stratification using MDS generated covariates. A significant improvement in the initially observed deviation from normal distribution of all 
p values as a result of correcting for relatedness and population stratification effects is observed. 



statistically significant markers identified in our GWA anal- 
ysis. An additional marker (BICF2P1279787) was identified 
in association with this haploblock (r2=0.87) with a reduced p 
value (P =0.00134). The 0.25 Mb locus on Chrl4 contains 

v genome ' 

20 SNPs and 2 CanFam2.0 annotated protein-coding genes, 
corresponding to the human genes COL1A2 and PEG10. 

A second statistically significant association was identi- 
fied (P =0.00049, OR=3.93, CI =1.78-8.66) for marker 

v genome 3 95% J 

BICF2P893476 on chromosome 24 (Figure 3). The marker is 
part of a two-SNP locus (Chr24: 43,091,222-43,595,979bp) 
spanning 0.5 Mb. RAB22A, a member of the rat sarcoma 
(RAS) oncogene family is located within the locus brack- 
eted by the two markers, which display complete LD (D — 1, 
r 2 =l). A notable difference in allele and genotype frequency 
distribution for BICF2P893476 was observed among cases 



and controls, at 37.8% "G" in cases versus 13.4% "G" in 
controls (Table 2). The allele and genotype distribution of 
BICF2P893476 is representative of all SNPs within the locus 
identified. 

In addition, five isolated SNPs displaying statistically 
significant P values were identified (Table 3). All vari- 

° genome v J 

ants identified are located in noncoding regions and were not 
further investigated. 

Tissue processing and histological analyses: Immunohis- 
tochemical expression of the two statistically significant 
PACG-associated genes, COL1A2 and RAB22A, was assessed 
in the dog eye. Staining of COL1A2 and RAB22A in sagittal 
sections of paraffin-embedded eyes revealed localization of 
both proteins in various substructures when compared to a 
negative control where the primary antibody was omitted 
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Figure 3. Manhattan plot from the 
genome-wide association (GWA) 
analysis of primary angle-closure 
glaucoma (PACG) in the basset 
hound (BH) with pedigree-based, 
multidimensional scaling analysis 
(MDS)-clustered covariates shows 
a significant association on chromo- 
some 14. The red line indicates all 
p values exceeding the permutation 
adjusted significance threshold of 
1(T 4 . 
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Table 3. Additional SNPs identified in GWA analysis. 
Chromosome SNP Position (bp) P 

v 1 ' genome 

26 BICF2G630805530 13,936,505 4.5x10-" 

20 BICF2P1094598 9,595,282 6.1x10-" 

32 BICF2P536453 21,851,662 6.2x10-" 

13 BICF2S23423784 58,788,894 6.5x10-" 

10 BICF2P528613 62,472,085 7.4x10-" 

Five single SNPs displaying statistical significance were identified in our case-control GWA analysis of PACG in the BH. 



(Figure 4). Prominent COL1A2 localization was noted in the 
ciliary body (CB) as well as the apical epithelial cell surface 
of the ciliary processes (Figure 4A). RAB22A similarly local- 
ized to the ciliary processes, but revealed faint staining in 
the CB, indicating reduced localization here. Both proteins 
displayed minimal localization in the CC and TM (Figure 
4B). Localization of COL1A2 was further noted in close 
proximity to the pigmented epithelium at the apical surface 
of the iris, whereas RAB22A staining was mostly observed 
in the stromal region (Figure 4C). In the cornea, COL1A2 
was found predominantly located in the stromal regions, 
whereas, RAB22A staining revealed strong localization in 
the epithelial layer at the apical surface (Figure 4D). Both 
proteins revealed localization throughout the retina. COL1A2 
showed a slight increase in staining in the outer nuclear layer 
and retinal ganglion cell layer. RAB22A revealed a strong 
staining signal in the outer segment and outer plexiform 
layer (Figure 4E). COL1A2 and RAB22A were also found 
to localize in the optic nerve when compared to the nega- 
tive control (Figure 4F). No obvious differences were noted 
in COL1A2 or RAB22A localization or intensity in tissue 
derived from affected versus unaffected dogs. 

To assess collagen localization and discriminate between 
type I and type III collagens in ocular tissue, Picrosirius 
red staining of paraffin-embedded BH eye sections was 
conducted (Figure 5). Type I collagen indicated by orange, 
pink, and red hues was the predominant collagen type 
observed in all structures of the eye. When compared to a 
negative control stained with Mayer's Hematoxylin only, 
abundant localization of relatively thick, red-colored, type 
I collagen fibers was noted in the CB and the epithelial 
apical surface of the ciliary processes (Figure 5A). Similarly, 
strong type I collagen staining, appearing in red, was noted 
throughout the iris, cornea, and optic nerve (Figure 5B,C,D). 
On the other hand, reduced abundance in collagen deposition 
and fiber thickness were noted in the CC and TM, as indi- 
cated by the apparent orange-pink-colored fibers (Figure 5E). 
This finding validates the observed reduction in COL1A2 
staining and localization in the TM shown above. Despite the 



absence of obvious differences in the localization or intensity 
of collagen staining in affected versus unaffected derived 
tissue, an overall widespread localization of collagen was 
noted in all eye structures, thereby confirming its prominent 
expression in the eye. 

DISCUSSION 

Despite recent discoveries, the underlying genetic and 
environmental factors that contribute to the development 
and progression of glaucoma have not been fully elucidated. 
The condition may result in debilitating outcomes, including 
complete loss of vision in the absence of effective treatment. 
The identification of genetic variants that determine risk for 
developing glaucoma can not only be used to genetically 
identify dogs at low disease risk for breeding purposes, but 
may also provide insight into the pathophysiology of PACG. 

GWA analysis of clinically confirmed BHs with PACG 
cases versus controls has led to the identification of two asso- 
ciated regions, which supports the role of a genetic component 
in modifying risk for PACG in the BH. The highest P 

Jo ° genome 

value was achieved for a susceptibility locus on chromosome 
14 containing the genes COL1A2 and PEG10. A second locus 
containing the gene RAB22A was identified on chromosome 
24. SNPs marking both loci achieved P values in the 

° genome 

order of 10 4 , which indicates significant despite the modest 
sample size used. The observed pattern of haplotype distribu- 
tion and allele frequency of markers within the loci identified 
does not support an additive effect of allele dosage in cases 
versus controls at multiple loci (Table 2 and Table 3). Based 
on these results, it appears that variation in multiple associ- 
ated regions determines risk for canine PACG rather than a 
single, major locus. These observations highlight the complex 
nature of the genetic and environmental effects underlying 
PACG inheritance. 

One shortcoming of this study is that in an effort to iden- 
tify a sufficient number of clinically well-characterized dogs, 
animals from multiple PACG pedigrees were used. Addition- 
ally, cryptic relatedness is particularly common for dogs; a 
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Figure 4. Immunohistochemistry staining for COL1A2 (left), RAB22A (middle), and negative control (primary antibody omitted; right) 
in sagittal sections of paraffin- embedded basset hound (BH) eyes derived from affected and unaffected animals. Variable intensities of 
COL1A2 and RAB22A staining were noted in the A: ciliary body (CB) and processes, B: trabecular meshwork (TM) and ciliary cleft 
(CC), C: stroma and pigmentary epithelium (PE) layers of the iris, D: epithelium (EP), stroma and endothelium (ED) layers of the cornea, 
E: retina (OS=outer segment, ONL=outer nuclear layer, OPL=outer plexiform layer, INL=inner nuclear layer, IPL=inner plexiform layer, 
GCL=ganglion cell layer; E), and F: optic nerve axons (ONA). Arrows indicate areas where staining was prominently observed. 
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Picro-sirius Red Negative Control 




Figure 5. Picrosirius red (left) and negative control-Mayer's hematoxylin only (right) stained sagittal sections of paraffin-embedded basset 
hound (BH) eyes. Images were viewed using polarizing light microscopy. Collagen staining of fibers appears in hues of orange, pink, and 
red. The observed color intensity is indicative of collagen fiber thickness and abundance in tissue. Collagen fiber staining was noted in A: the 
ciliary body (CB) and processes, B: stroma and pigmentary epithelium (PE) layers of the iris, C: epithelium (EP), stroma and endothelium 
(ED) layers of the cornea, D: optic nerve (ON), E: ciliary cleft (CC) and trabecular meshwork (TM). Arrows indicate areas where staining 
was prominently observed, including the CB, iris, lens, and ON. 
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certain amount of relatedness is almost always detected, 
particularly in purebred dogs [22-24]. Sample independence 
is one of the main assumptions of most GWA test statistics 
[25,26], and thus the degree of relatedness among animals 
is a common concern in GWAS. However, a large number 
of animal studies were able to come to valid conclusions 
despite the discovery of cryptic relatedness among their study 
samples following the investigation of population stratifica- 
tion [22,23,27]. Herein, sample relatedness was addressed by 
selecting cases and controls that were as distantly related to 
each other as possible and by using statistical approaches to 
correct for confounding relatedness. Our results reveal the 
effectiveness of the statistical methods used here to resolve 
the confounding effects of cryptic relatedness initially 
observed among our samples. 

Our findings also demonstrate that significant pheno- 
type-genotype associations can be achieved with a modest 
sample size when using a closed and purebred population, 
such as the BH population used here. When compared to 
human populations, studies of the canine genome have shown 
that genetic homogeneity is much greater within individual 
dog breeds (94.6%) than within distinct human populations 
(72.5%) [28]. In some breeds, genetic variation has even been 
additionally reduced by extensive inbreeding and bottleneck 
events. The vast reduction in haplotype diversity in the dog 
genome facilitates mapping the genetic basis of phenotypic 
variation using GWAS [29,30]. Genetic associations can 
therefore be identified more efficiently using fewer markers 
and samples in canine versus human genetic studies [31]. 
With that in mind, the small sample size used in this study 
may be considered too modest to provide sufficient statistical 
power for detecting small genetic effects. This is particularly 
true if the risk variant is a rare allele with a small OR [32]. 
Human-based studies have demonstrated that when the causal 
SNP is rare (minor allele frequency [MAF]<10%), a sample 
size of 3,000 may still be underpowered to detect small 
genetic effects [33]. These findings may not be applicable to 
canine -based genetic studies considering the architecture of 
the canine genome, which could provide greater statistical 
power for the identification of risk variants with small effects 
using a smaller sample size relative to human-based studies 
[28]. 

It has been shown that many loci harboring common risk 
alleles for complex diseases will have effect sizes in the range 
of 1.1-1.2, which cannot be detected using the conventional 
study designs and sample sizes currently adopted, but rather 
using empirical analyses and significantly large sample sizes 
[34,35]. As such, our inability to detect small genetic effects 
with the sample size used here is expected considering that 
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larger GWAS remain statistically underpowered to detect 
common variants with small effects. Our results reveal the 
strong association of two loci, containing promising candi- 
date genes, to PACG. These findings support our expectation 
of detecting variants with large effects in association with a 
complex disease when using a small sample size. 

Several research studies have provided convincing 
evidence of the potential role of collagen in glaucoma, making 
COL1A2 a highly promising PACG candidate gene. Of partic- 
ular interest is that a variant in the collagen gene COL11A1 
has been recently identified in a GWAS of human PACG [8]. 
Central corneal thickness (CCT) is a known and proven risk 
factor of glaucoma [36]. In a GWA study investigating CCT 
and brittle cornea syndrome in association with glaucoma, 
significant associations in the collagen genes COL8A2 and 
COL5A1 were identified [37]. These findings suggest a 
possible involvement of collagen genes in the pathogenesis 
of PACG. Moreover, analysis of AH collected from patients 
with POAG revealed increased fibroblast proliferation and 
collagen synthesis in POAG versus unaffected derived AH 
[38]. The increased rate of fibroblast activity is speculated to 
contribute to excess deposition of collagen and subsequent 
loss of the TM cells during the development of glaucoma 
[38]. In addition, this locus contains a second gene — PEG10. 
Little is known about this gene, but studies in the mouse 
demonstrate that it is an imprinted, retrotransposon-derived 
gene. PeglO knockout mice showed early embryonic lethality, 
which indicates its critical role in mouse development [39]. 
Ras-related Protein RAB22A, a member of the RAB family 
of small GTPases was identified within the second PACG 
associated locus in the BH. The protein encoded by this gene 
has been shown to interact with early-endosomal antigen 1, 
and may be involved in the trafficking of molecules between 
endosomal compartments [40]. Neither PEG10 nor RAB22A 
has been discussed in relation to the pathology of glaucoma 
in the past. 

In conclusion, we have identified two loci associated 
with PACG in the BH. Our findings indicate the possible 
involvement of a collagen-related mechanism in the develop- 
ment of glaucoma in the BH. This finding is in accordance 
with emerging data that support the role of collagen in human 
PACG. Additional functional studies will be required to 
investigate the role of COL1A2 and possibly other collagen- 
encoding genes in PACG. This GWA analysis of genetic 
associations has demonstrated the complexity of PACG 
by implicating several genetic loci in association with the 
disease. Sequencing of candidate genes within the identified 
disease-associated loci will possibly lead to the identifica- 
tion of PACG-causing mutations in the BH. In addition to 
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identifying the genetic variants underlying the disease, 
understanding the environmental factors that potentially 
modify risk to glaucoma is crucial for characterizing the 
pathophysiology of PACG. Ultimately, we anticipate the utili- 
zation of our findings in estimating disease risk in susceptible 
BHs and possibly other dog breeds susceptible to PACG. This 
finding may contribute to the reduction of disease incidence 
in the BH by selectively breeding animals following genetic 
testing for risk variants. With further research, findings in 
this study may shed light on important molecular mecha- 
nisms that contribute to PACG in the canine and potentially 
human patients. 
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